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Abstract 

The Skyrme model, an effective low energy theory rooted in large N c QCD, has been applied 
to the study of dense matter. Matter is described by various crystal structures of skyrmions. 
When this system is heated, the dominating thermal degrees of freedom are the fluctuating 
pions. Taking these mechanisms jointly produces a description of the chiral phase transition 
leading to the conventional phase diagram with critical temperatures and densities in agreement 
with expected values. 
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An important issue today is to understand the properties of hadronic matter under extreme 
conditions, e.g., at high temperature as in relativistic heavy-ion physics and/or at high density 
as in compact stars. The phase diagram of hadronic matter turns out richer than what has 
been predicted by perturbative Quantum Chromodynamics (QCD) pQ. Lattice QCD, the main 
computational tool accessible to highly nonperturbative QCD, has provided much information 
on the the finite temperature transition, such as the value of the critical temperature, the type 
of equation of state, etc[2]. However, due to a notorious 'sign problem', lattice QCD could not 
be applied to study dense matter. Only in the last few years, it has become possible to simulate 
QCD with small baryon density [3]. Chiral symmetry is a flavor symmetry of QCD which plays 
an essential role in hadronic physics. At low temperatures and densities it is spontaneously 
broken leading to the existence of the pion. Lattice studies seem to imply that chiral symmetry 
is restored in the high temperature and/or high baryon density phases and that it may go hand- 
in-hand with the confinement /deconfinement transition. The quark condensate (qq) of QCD is 
an order parameter of this symmetry and decreases to zero when the symmetry is restored. 

The Skyrme model, an effective low energy theory rooted in large N c QCD, has been applied 
to the dense matter studies [U El El El El [TO] . The model does not have explicit quark and 
gluon degrees of freedom, and therefore one can not investigate the confinement /deconfinement 
transition directly, but we may study the chiral symmetry restoration transition which occurs 
close by. The main ingredient associated with chiral symmetry is the pion, the Goldstone 
boson associated with the spontaneously broken phase. The various patterns in which the 
symmetry is realized in QCD will be directly reflected in the in-medium properties of the pion 
and consequently in the properties of the skyrmions made of it. 

The classical nature of skyrmions enables us to construct the dense system quite conveniently 
by putting more and more skyrmions into a given volume. Then, skyrmions shape and arrange 
themselves to minimize the energy of the system. The ground state configuration of skyrmion 
matter are crystals. At low density it is made of well-localized single skyrmions At a 

critical density, the system undergoes a structural phase transition to a new kind of crystal. It 
is made of 'half-skyrmions' which are still well-localized but carry only half winding number. In 
the half-skyrmion phase, the system develops an additional symmetry which leads to a vanishing 
average value for (a) [12]. In the studies of the late 80's [13], the vanishing of (a) was interpreted 
as chiral symmetry restoration by assuming that (a) is related to the QCD order parameter (qq). 
However, in Ref. [1], it was shown that the vanishing of (<r) cannot be an indication of a genuine 
chiral symmetry restoration, because the decay constant of the pion fluctuating in such a half- 
skyrmion matter does not vanish. This was interpreted as a signal of the appearance of a 
pseudo-gap phase similar to what happens in high T c superconductors |14j . 

The puzzle was solved in Ref. [7] by incorporating a suitable degree of freedom, the dilaton 
field X) associated to the scale anomaly of QCD. The dilaton field takes over the role of the order 
parameter for chiral symmetry restoration. As the density of skyrmion matter increases, both 
(a) and (x) vanish (not necessarily at the same critical density). The effective decay constant 
of the pion fluctuation vanishes only when (%) becomes zero. It is thus the dilaton field which 
provides the mechanism for chiral symmetry restoration. 

Contrary to lattice QCD, there are few studies on the temperature dependence of skyrmion 
matter. Skyrmion matter has been heated up to melt the crystal into a liquid to investigate 
the crystal-liquid phase transition [151 [TB] a phenomenon which is irrelevant to the restoration 
of chiral symmetry. In here, we study skyrmion matter at finite density and temperature. We 
obtain the phase diagram describing the realization of the chiral symmetry. 

The basic strategy of our approach is as follows. We describe dense skyrmion matter as 
a certain crystal structure. Thereafter, we let the temperature rise. The dominating thermal 
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degrees of freedom in this process are the fluctuating pions which contribute a thermal energy 
proportional to T 4 . We neglect the thermal fluctuations of the heavy skymions (translations, 
rotations and vibrations) which are proportional to T. We show that this strategy leads to a 
familiar phase diagram with critical temperature and density in agreement with the expected 
values. 

The simplest model Lagrangian containing the inevitable ingredients to make the calculation 
respect the flavor properties of QCD reads [18] 

£ = X 2 £ a + £ sk -V( X ), (1) 
with ^ 

6/e sk 

V(x) = ^{x 4 (lnx -|)+|}. (fffc) 

In this equation, U = exp(zr • tt//-k)(& SU(2)) is a nonlinear realization of the pion fields and 
X represents the dilaton field. The model Lagrangian contains four parameters, /„■, e sk , m x and 
f x , i.e. the pion decay constant, the Skyrme parameter, the dilaton decay constant and the 
dilaton mass, respectively. We choose in what follows f n = 93 MeV, e s k = 4.75, m x = 720 
MeV and f x = 210 MeV (a discussion on the values of these parameters can be found in 
Ref. |7j). For simplicity we neglect the pion mass and take the x field as a constant in space 
and time. Still, the 'constant' dilaton field controls the chiral phase transition through the 
potential energy term. In the vacuum, this term has its minimum at x = 1. At finite density 
and temperature, skyrmions and thermal pions contribute to the term proportional to x 2 an d 
the system undergoes a dynamical phase transition when this term is larger than the potential 
energy. While the specific values for the critical temperature and density depend strongly on the 
chosen parameters, the phase transition scenario remains qualitatively the same for all reasonable 
values. 

Skyrmion matter is described by the energy of a single skyrmion: 

£(p) = x 2 £* + 4k + V( X )V, (2) 

where 

S a (p) = [ d 3 r tr(diU diU^, m 
Jv 

£sk(p)= / d 3 r tvlU^diUo^djUo] 2 , d2b) 

p is the baryon number density of the system and V is the volume occupied by a single skyrmion; 
i.e., p = 1/V. Uq{t) is the static skyrmion configuration that minimizes the system energy £ 
for a given p. Note that both the potential energy terms, which contains an explicit \/p{= V) 
factor, and £ CT)S k, through Uq{t), depend on the baryon number density. 

What happens if we heat up the system? Naively, as the temperature increases, the kinetic 
energy of the skyrmions increases and the skyrmion crystal begins to melt. The kinetic energy 
associated with the translations, vibrations and rotations of the skyrmions is proportional to T. 
This mechanism leads to a soliddiquid-gas phase transition of the skyrmion system. However, 
we are not interested in this transition but in the chiral symmetry restoration transition. A new 
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mechanism must be incorporated to describe it and it happens to be the thermal excitation of 
the pions in the medium. This phenomenon is proportional to T and therefore dominates the 
absorption of heat. We only consider the latter in our discussion. 
The pressure of the non-interacting pions is given by [19] 

P = — T 4 , (3) 
30 ' y ' 

where we have taken into account the contributions from three pions, 7r + ,7r°,7r~. This term 
contributes to the energy per single skyrmion volume as 3PVx 2 - The factor 3 comes from the 
fact that our pions are massless. The kinetic energy of the pions arises from C a , and therefore 
scaling implies that it should carry a factor x 2 ■ After including the thermal pions, the energy 
lodged around a single skyrmion can be expressed as 

£(P, T) = \S a (p) + ^T 4 yj x 2 + £ sk (p) + V( X )V. (4) 

In order to complete our calculation we have to determine the values of Uo(r) and x that 
minimize £(p,T) for a given density and temperature. 

Chiral restoration will occur when x vanishes. How does this take place in hot and dense 
matter? In eq.Q, the potential energy of the dilaton is finite and its magnitude decreases as 
the density p increases (V decreases). The potential energy has a potential hill at x = and a 
valley at % = 1- Beyond x = 1) h increases monotonically. Besides the potential term there is 
another term proportional to x 2 ■ As the density increases the S a contribution becomes relatively 
important. The thermal pion contribution is proportional to T 4 /p. Thus, the term proportional 
to x 2 increases to overcome the potential hill. 

Let us look into this phenomenon in detail. For a given p and T, explicit dependence of 
eq.(JI|) on x can be presented as 

e = ax 2 + b + c( X \\nx-\) + \)- (5) 

One may easily match the coefficients a, b and c with eq.([4]). S has two local minima at x = 
(chirality symmetric phase) and at x = Xo(t^ 0) (chirality broken phase) satisfying 

a + 2c X 2 lnxo = 0. (6) 
The energy difference AS = S(x = 0) — S(x = Xo) is obtained as 

AS = \xl{2a-cxl). (7) 
The phase transition happens when AS = 0; that is, 

2a - cxl = 0. (8) 
Combining eqs.([6]) and flSJ), we get xo at the critical point as 

X c o = e" 1/4 . (9) 

Note that the phase transition happens suddenly from a non- vanishing x = Xo value to x = 0- 
It is therefore a first order phase transition. However, this result is just a peculiarity of the used 
dilaton potential V(x)- If we had taken V{x) ~ (x 2 ~ I) 2 ) the phase transition would have been 
second order. 
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Using the above equations at the critical point, we get 



tt 2 „, flm\ 



which leads to 



^ + io^ = &f- (10> 



For p = (zero density), our estimate for the critical temperature is 

/in f 2 m 2 \ 1/4 

where we have substituted f x = 210 MeV, m x = 720 MeV. It is remarkable that this naive toy 
model leads to T c ~ 200 MeV, which is quite close to that obtained by lattice QCD [2] and in 
agreement with the data [20J. 

To complete our study we compute E a as a function of p by minimizing numerically £(p, T) 
with respect to Uq{t) and x- An useful approximation may be obtained by fixing \ = 1 and 
minimizing £ a (p) + £ s k(p) with respect to Uo(r). Since \ varies only in the range e _1//4 ~ 0.8 < 
X < 1, before the system undergoes the phase transition, the approximation is quite reasonable. 
These minimum energy configurations can be found by various methods [21], and they turn out 
to be single skyrmion FCC (Face Centered Cubic) crystals at low density and half-skyrmion CC 
(Centered Cubic) crystals at higher density. 

The dependence of £ a and £ s k on f n and e s k can be separated out explicitly as 

£<7,sk(/°) = (U/eak)i cr ,ak(p), (13) 
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Figure 2: The chiral phase transition. The solid line shows the exact calculation, while the gray 
lines two approximate estimates. 



where p = p/ie^f^) 3 and £ a , s k are dimensionless quantities which do not depend on e s k and 
/„-. Shown in Fig. 1 are the numerical values on £ a and In order to show their dominant 
asymptotic behaviors, we present various values of £ as a function of l/{p) 1 ^ 3 . For high densities 
£(j is linear in l/(/5) 1 / 3 , while for low densities it stays constant, implying that at a low density 
the interaction between skyrmions is negligible. We parameterize £ a as 



£c 



10/2/^ p> Po 

36/Tr/esk, p < po, 



(14) 



where p = (e sk f n /3.6) 3 . 

Using eq. (|14p for £ a and returning to eq. (jl(jp we obtain the critical density for chiral 
symmetry restoration at zero temperature as 

/ f 2 m 2 1 \ 3/2 

p c (T = 0) = "TTT^TTTT? ~ °- 37 fm ^- ( 15 ) 



k 8eV2 10/2 J 

Since po = 0.24 fm~ 3 < p c (T = 0) we have been consistent when using the high density formula 
for £fj 

The resulting critical density p c (T = 0) ~ 0.37 fm~ 3 is only twice normal nuclear matter 
density and it is low with respect to the expected values. This result is not problematic since 
p c (T = 0) ~ (fx m x/f^) 3 an d T£ =0 ~ (/x 771 x) 1//2 an< ^ therefore the description of the phase 
transition does not depend on the specific values of the parameters. Moreover, small changes in 
them may lead to larger values for the critical density. 

For a finite density smaller than p c (T = 0), we obtain the corresponding critical temperature 
by substituting the asymptotic formulas (fT4"l) for £ a , 
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T c = 
T c = 



T c (p = 0) (1 - 3.09(/m 3 ) p c ) l ' A for p < p 
T c (p = 0) (1 - 1.92(/m 2 ) p2/3)V4 for p>pQ 



(16) 
(17) 



The results of these two curves are shown in Fig. 2 by gray lines. 

The exact calculation (black dots connected by black line in Fig. 2) has been obtained 
numerically by minimizing the energy eq. ([!]). The resulting phase diagram has the same shape 
but the temperatures and densities are generally smaller than in the approximate estimates 
shown in Fig. 2. 

There is a long history of success of the Skyrme model in describing the chiral symmetry 
restoration phase transition for dense hadronic matter [U O El El El El ID] • The model might 
also be a viable description of the new phenomenology that is being proposed for dense/hot 
matter |22j . The time had come for an analysis of its behavior with temperature. We have 
presented a first description of the chiral restoration phase transition in the temperature-density 
plane, whose main ingredient is that the dominant mechanism is the absorption of heat by the 
fluctuating pions in the background of crystal skyrmion matter. This description leads to a phase 
transition whose dynamical structure is parameter independent and whose shape resembles much 
the conventional confinement /deconfiment phase transition. Moreover, for parameter values 
close the conventional ones, we obtain the expected critical temperatures and densities. Further 
investigation in these matters is required since it is becoming apparent that the path to the 
phase transition is not without new physics and that this 'bottom up' approach might be useful 
to obtain interesting new physics [23J. 
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